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\ Portnoi and Galbraith recently proposed a beautiful and intriguing relationship defining the crit- 

■ ical screening lengths associated with the apparition of new bound states for the two-dimensional 

' statically screened Coulomb potential. Not only does semiclassical quantum theory show that this 

relationship is unfortunately not strictly exact, it has also proved helpful in the search for a poten- 

= . -e^ 

, tial which exactly verifies Portnoi and Galbraith's formula, namely — -. — ttt. The analytical 

— » r{l+r/rsY 

eigenfunctions of the corresponding Schrodinger equation at zero energy, when localized, lead to 
approximate upper bounds for critical screening lengths of the two-dimensional statically screened 
Coulomb potential. 
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I. INTRODUCTION 

The two-dimensional statically screened Coulomb potential Vs (r) plays a central role in the physics of semiconductor 
heterostructures near the band gap, especially with respect to their linear and nonlinear optical properties when the 
influence of excitons is properly taken into account Although already considered three decades ago it 

' has nonetheless received a lot of recent attention and been analyzed using different approaches such as the WKB 
I approximation J0| , perturbation theory |^ , variational calculation [^[^ or numerical resolution based on the variable- 
O ' phase method |pPpT| . 

O , Using a variable-phase approach ^ to count the bound states allowed by Vs{r) — only a finite number of them 
exist — and making a connection with the Levinson's theorem in two dimensions, Portnoi and Galbraith recently 
proposed a beautiful and intriguing relationship between the critical screening lengths (scaled to the exciton Bohr 
^ radius) and the number of bound states for a given value of the angular momentum (Z > 0): 

^ ■ . . . 

(the original notation used v = iV + 1, whereas we focus on the number N of nodes of the wavefunctions). They 
T-H proceeded to derive a formula giving the number of bound states as a function of Ts , which markedly differs from a 
WKB estimate |^ and from the Bargmann bound condition reformulated for the two-dimensional case [ pl| . 

In this Communication, we show that semiclassical quantum theory simply demonstrates that Eq. (|l|) is not strictly 
exact for Z = states, and that the previous estimate using the WKB approximation is incorrect. We then explain 
y ' how this semiclassical approach has also proved useful in the search for a potential which obeys Eq. ([|), for which 
\ \ the eigenfunctions at zero energy are analytical. Finally, we show why, for the states defined by iV = and I > 2, the 
(— ( true critical screening lengths of Vs{r) are nearly given by Eq. (|l|). 

o . 
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II. SEMICLASSICAL QUANTUM THEORY AND THE TWO-DIMENSIONAL SCREENED COULOMB 

POTENTIAL 

The two-dimensional screened Coulomb potential Vg has a very simple form when expressed in wavevector space 

VM = -^, (<ls^-) (2) 
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where we have adopted the same units as in Ref. |^ (the length scale is the effective Bohr radius gq and the energy 
scale is the three-dimensional exciton energy R). In real space, it reads [0j|J§l 

Vsir) = -^[l-^q,r[lio{qsr)~No{qsr)]} , (3) 

where Hg and Nq are the Struve and Neumann functions, respectively. The term between braces in Eq. goes 
down to like /r^ as r goes to infinity |^ , which is a much reduced decrease with respect to the exponential decay 
occurring in three dimensions. 

Let us now turn to the semiclassical quantum theory approximation for a potential V{r). Using conformal mapping, 
Yi and collaborators have shown [|l2| that the two-dimensional formulation of energy quantization may be written as 



[N > 0) (4) 



where Vcs{r) = V{r) + P/r'^, while ri and r2 are the classical turning points. They have proved that Eq. (|4|) gives 
exact results in several cases, and excellent ones for the energy levels of impurity states in an arbitrary external 
magnetic field. Since we are looking for the critical values of at which new bound states appear at zero energy, we 
can set i? = in Eq. (0) so that, for the I — states of V^, 



(5) 



In the following we call (rg)^ the value obtained by using the semiclassical approximation in contrast to the exact 
one. Eq. (||) gives 



irs). = -,[^) {N+-^] (6) 



I=J^ ^l--u^Ilo{u^)-No{u^)]du. (7) 

As is usual with semiclassical expressions, we expect Eq. ^ to be accurate when A'' ^ 1. If Eq. (P were true, we 
would find 21 = TT, so as to get the correct leading dependence of {rs)c in N'^- ^ has been numerically evaluated 
by first calculating the integral between and 10, and then bracketing the integral from 10 to -|-oo by using the 
asymptotic expansion of Hg — Nq. The same procedure, repeated for an intermediate bound of 20, gave the same 
result, namely 

2X^ 3.14057414. (8) 

It definitely differs from tt, even though the relative error is only about 3.24 x 10^**. Consequently, for ^ = states 
and large A'''s, (rg)^ should be different from Eq. (|l|). 



Let us now discuss previously published results in the light of Eqs. (|6|) and (g). Portnoi and Galbraith JLO 11 
compared Eq. (Q) with the WKB prediction of Reyes and del Castillo-Mussot for I — states and observed 
a disagreement by a factor 2.5. By contrast, Eq. (||) satisfactorily removes the discrepancy. Why then do our 
semiclassical approximation and the WKB result differ by such a large amount ? The explanation probably comes 
from the tricky use of WKB wavefunctions and asymptotic expansions in Eqs. (7)-(12) of Ref. In particular, 
their Eq. (12) contains w"^/^ vT^^^lw dw, which diverges. Admittedly, this integral can be formally written as 
5(3/2, —1/6), whose true value is -6.72 instead of the given -2.81. This oversight seems to have been left uncorrected 
until the final result in their Eq. (18), casting therefore strong doubt about its validity. 

At this point, we have merely shown that Portnoi and Galbraith's formula cannot be exact for I = states, since 
it does not verify the large N limit. Strictly speaking, we cannot use this argument to reject their hypothesis that 
(rs)c =1, 3, 6, 10, ... for I — and A*" = 1, 2, 3, 4 ... or more generally, for low-energy states. However, because of 
the near-coincidence of 2T with tt, one may have second thoughts about the accuracy of their (rg)^: the statically 
screened Coulomb potential is very long range because of its l/r'^ decrease at infinity |jl^ . Besides, a question remains 
unanswered: is it possible to understand why the critical screening lengths are so close to the integers of Eq. (ml 

2 r 

For this reason, we investigated what happens if Vg is replaced by — /( — ), when / ^ 1 for r ^ and f(x) 

r rg 

monotonically decreases as l/x^ for x oo. 
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III. SEMICLASSICAL APPROXIMATION FOR A FAMILY OF SCREENED COULOMB POTENTIALS 



2 r 1-1 

For an attractive potential — /( — ), Eq. (|]) written at zero energy gives 

r Ts 



(^s).-^-^ ^ (^iV+^l . (9) 





Let us call A the prefactor of {N + 1/2)^. We see on Table | that even though the asymptotic behaviors are the same 
for all the potentials (— 2/r and —2 r^/r^ for r going to and oo, respectively), A may vary between at least 1/3 and 
2/3. The criterion of apparition of bound states is thus extremely sensitive to the exact shape of the potential over 
its whole range of variation. Obviously, for 

yair) = -l . (10) 
r f r 
+ — 

we have A = 1/2. Because of its simple form, one can further investigate this potential for I ^ 0. The calculation is 
performed without difhculty and yields 

Equation (^Tj) strikingly resembles Eq. (]l|), apart from a difference of 1/8. Since simple results derived from semi- 
classical quantum theory are often closely related to solvable problems, we have searched for the exact analytical 
expression of the critical screening lengths for Va given in Eq. ([To|). 

IV. EXACT CRITICAL SCREENING LENGTHS FOR Va 

For the potential Va, the Schrodinger equation verified by the radial wavefunctions ipo{r) of zero energy reads 

< (0 + - ^0 W + ( - - 4 ) Mr) = 0. (12) 

r \r(rs+ ry j 

Because of the particular form of T4, one may first consider rational fractions of r for -00. Starting with = 
(for wavefunctions having no node), one can check that r'- /{r + Tj)^' obeys Eq. (|l2|) provided that Ts — I {21 + 1). 
Proceeding further with N ~ 1,2,..., it is not difficult to see that the structure of the solutions is 



2r 



Mr) « -IVT2T -^:-^^-2^;2/ + l;- , (13) 

/ j {N+2i)(N+2i+i) Y^ \ [N + 21) [N + 21 + I) ' 

where 2^1 is the hypergeometric function, which here reduces to a polynomial of degree N . Quite remarkably, the 
critical screening lengths are exactly given by Portnoi and Galbraith's formula, Eq. (|l]). This also shows a posteriori 
that semiclassical quantum theory can be trusted to provide the correct result when N is large. 

A more direct way to solve Eq. dl2) is to set ipoir) — r' x( ); x(") then satisfies the hypergeometric equation. 

Eq. (|l|) is obtained through the condition that the radial wavefunction, which is finite at the origin, does not diverge 
at infinity: one of the first two arguments of 2-F1 must then be a nonpositive integer. This leads to 

Mr) oc r' ,F,i^^ N ^ 21 N + 21 + 1,21 + 1,—^^^^^^^^ , (14) 

which is equivalent to Eq. ( |l3| ) p[ . These wavefunctions are currently considered for the determination of upper 
bounds for the critical screening lengths of Vs , when I > 1 . 
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V. UPPER BOUNDS FOR CRITICAL SCREENING LENGTHS OF K 



The above wavefunctions have an interesting property for I >2, which is obvious from Eq. (]l3|): they are locahzed 
(square integrable). Since Va and Vg have the same asymptotic behaviors and seem to have nearly identical (rg)^, we 
may be tempted to use perturbation theory to get an estimate of the critical screening lengths of V^. In the foUowing, 
we will restrict the discussion to the {N = 0, / > 2) states. 

Starting from the exact solution ipQ of zero energy given in Eq. (|l3|), the first-order correction to the energy is 
obtained by calculating the expectation value of Ay = Vg — Va on this state, namely 

^<=o,>2 = (^ol VAr) + - IV'o) . (15) 

(r + /(2/ + l)) 



Assuming that (r^)^ — I {21 + 1) for Vs leads to (after normalization of tAq) 



^N=0,l>2 



;(2/ + l)(2Z + l)!(2Z-3)! 



X 







The above integral is much easier to evaluate numerically than T, because m! (1 + u) ™ ^ is the Laplace transform 
of f " e"* and because Q 

1 , /(i + vm^)(i + \/rTi^)x 

rfu - [Ho(m) - iVo(u)] e "* = ^- In ^ ^-^ . (17) 







One finds that Ai^j^^^g is equal to 1.29 x 10""*, 1.22 x 10"*, 9.87 x 10"^ 7.83 x 10~^ 6.27 x 10~^ 5.10 x 10"^ 
4.22 X 10"^ for I = 2, ...,8, respectively. These small positive values seem to indicate that for the lowest-lying states 
with Z > 2, the critical screening lengths of the two-dimensional statically screened Coulomb potential are probably 
larger than ^ (2 Z -t- 1), but only slightly so. 

A visual way to reach this conclusion is to consider the different functions appearing in the integral of Eq. (|T^). 
We see in Fig. ^ that r AV(r) goes to zero at the origin and r oo, as expected; however, its sign is not uniform. 
In order to get AiJ^^g ;~>2i o^ie must also consider the squared wavefunctions, |'0oP, the amplitude of which rapidly 
decreases with I. Furthermore, the maximum of IV'oP is always located at it = 1, i.e., r ^ I {2 1 + 1); obviously, positive 
and negative contributions to the integral nearly compensate each other. Since the first-order energy corrections are 
very small already, we can understand why the critical screening lengths for Vs must be very close to 1(21 + 1) for 
the [N ^0,l>2) states. 

For these states, we can actually use second-order perturbation theory to get an admittedly rough upper limit to 
the exact value of (r^)^ for Vs. The argument is the following: it is well known that second-order energy corrections to 
the state of lowest energy are always negative. This result can be extended here to the iV = states, since subspaces 
corresponding to different Vs remain uncoupled by the isotropic AV. If, instead of taking rs = l{2l + 1) as in 
Eq. (p^), one considers = l{2l + l)/ry/ with rji such that A_B^^q = 0' then A_B^^q < 0; which suggests the 
possible existence of another bound state and consequently (rg)^ < l{2l + l)/rji. As regards calculations, one merely 
has to replace aU the u's in the Vs part of Eq. (|l|) by rjiu. Finally, one gets (rs)!,'^^] < 10.106, (rs)[.'=3l < 21.273, 
(rs)[f^*' < 36.524. Although rather crude, these upper bounds have been verified (and improved) by a more accurate 
variational approach, which also showed that (rs)'^"^''"^^ is strictly less than 6 [Ol. 



VI. CONCLUSION 

In conclusion, we have shown that Eq. (|l|), proposed by Portnoi and Galbraith to describe the critical screening 
lengths for the two-dimensional statically screened Coulomb potential Vs, is unfortunately not strictly exact. This 
does not invalidate their variable-phase method approach, which is very useful: it is quite likely that the discrepancy 
between the exact values and the integers they found must be very small. Using analytical eigenfunctions at zero 
energy for a potential very similar to Vs , we have been able to explain why the critical screening lengths corresponding 
to the fundamental states for Z > 2 are indeed nearly given by Eq. ([^). 
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TABLE I. Values of A such that {rs)c = A {N H — ) in the semiclassical approximation, for different potentials /( — ). 
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TV L U 

FIG. 1 . Behavior of 1 u [Hn (u) — Nn (u)] — -. (full line) and of -. — -r , which are the squares of the wavefunctions 

used in Eq. ([l^). For clarity, the curves for i = 2, 3, 4 have been multiplied by 10, 100, and 1000, respectively. 
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